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Abstract 

Jeans mass calculated with different combination of parameters 
involved has shown interesting variation with remarkable shifting of 
position of peak value from x = m/T = 0.5 to 5.5. The standard 
deviation is 2.217. In particular, using the harmonic mean square ve- 
locity, shifts the peak Jeans mass to x = m/T ~ 2, which is remarkably 
less than previously reported value of 4.2. Different scales of neutrino 
structures including virialized moments have also been compared. 

Introduction 



The enigmatic neutrino has generated great deal of interest in cosmology for 
a long time. Regarding the number of neutrino species e.g., calculation of 
primordial nuclear abundances in cosmology had arrived at the number of 
neutrino species A*"!, < 4 [1] long before it was experimentally established at 
CERN|2]; the current value is iV^ = 2.994±0.013 [3]. Similar calculations set 
the baryon to photon ratio at ~ 4.7 x 10~^°, which in turn has determined 
the present baryon density in the Universe to be Q^/i^ = 0.0223 , where 
= Pb/Pc is the present baryon density in units of the critical density, and 
h is the present value of the Hubble parameter in units of 100 km/s/Mpc, 
that is expected to be /i ~ 0.72(3) |4]. Writing the present neutrino density 
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as Q^h'^ = mf.y/?)\, where is the average rest-mass of the three species 
neutrino in eV , it is easily seen that the 31 eV neutrino will close the Universe; 
this closure mass is referred to as the Cowsik-McClelland bound|5j, although 
in their paper they had used a smaller value of VL and four component neutrino 
to arrive at Se^. Tremain and Gunn [6] set the limit that a neutrino of 
rrijj < IMeV cannot reside in the halo of galaxy to contribute significantly 
to the dark matter. Bond et al [7] calculated the maximum Jeans mass for 
structures of neutrino to be M^m ~ 1-2 x 10^'^m~yMQ , where mev is the 
neutrino mass in eV. 

As recent results show that the mass difference between neutrino species 
is very slight, Am^ ~ 10^^ eV^, in this paper, we consider three neutrino 
species of similar mass. Comparing the baryon and neutrino densities, we see 
that even a 0.7 eV neutrino will have dominated over baryon by now. Rather 
than whether the neutrino is massive, the question at present is how mas- 
sive it is; again cosmology provides the most stringent limit "Y^rriv < 0.6 
[TT| [T2] . So the neutrino appears as a very important component of the Uni- 
verse. Consequently, it should have a strong bearing on structure formation. 
It decouples at very early times, and then evolves as a totally independent 
component that interacts only gravitationally. Indeed the filamentary struc- 
tures, sheets, walls and voids as exposed by various surveys point towards 
dissipationless collapse of clouds of particles like neutrinos into Zeldovich 
pancakes 1 13j at some stage of evolution. Although the cosmological neutrino 
has been studied extensively, we feel that our analysis gives some further in- 
sight on structure formation, particularly at what temperature the structures 
are more likely to form and spectrum of the size / mass of the structures at 
various temperatures. 

The momentum of a freely moving particle in the Friedmann-Robertson- 
Walker spacetime is redshifted by the expansion, i.e., the comoving momen- 
tum pa = y remains constant, where p is the momentum and a the scale 
factor. Using v = p/E for the velocity and the Einstein energy- momentum 
relation, a'^E'^ = p^a? + im?a? = + x^, we see that Eva remains constant 
during expansion, and as the number density scales as a~^, pva^ also re- 
mains constant where p is the density. As the light neutrinos decouple at the 
very high temperature of T ^ IMeV while still extremely relativistic (ER), 
they are essentially in free fall since then. So their number density is always 
distributed as 
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g fm\^ y"^ 



where g is the number of spin degeneracy (six, for the three u — V pairs), 
and we have used the fact that T ~ 1/a ; Planck units in which G = c = 
ks = h = 1 are used throughout this work. As ~ p >> m in the extreme 
relativistic (ER) regime, Eq.([T]) represents the Fermi-Dirac distribution. In 
the non-relativistic (NR) case however, E ~ /2m + m, and Eq. ([1]) is no 
more a Fermi-Dirac distribution. 

Integrating dn over y from to 00 gives the number density 
n = ^T{3)ri{3) (^) where F and 77 are the gamma and eta functions 
respectively |15|: also, 77(72) = (1 — 2^^"^) ({n). Thus we can write down the 

expectation value of any regular function as < f{y) >= r(3)',;(3) Io° ^V'^+ifiv) 
The width of the momentum distribution shown in Fig. [flis characterized by 



Umax = 3.131 where the distribution is maximum. Some other characteristic 
values that we will use are the mean ymean =< V >= rl^)r)^^) ~ Ic(t) ~ 3.151, 

the root-mean square yrm.s =< >^/'^= -y/^^^y^ = 3. 597, the harmonic 

mean yhm =< U^^ >^^= = 2.192, and the root harmonic mean square 



Vrhms =<y ^ > ^/^= = Y = 1-613 are also useful to describe 

the width of the distribution. 



The particle speed is v = p/E = y / y^ + x'^. In the ER case as y >> x. 



V{y) I - + giving Vmean ^ 1 " J (^;;f^) Vrhms] alsO, Vhni =< 

1/v >^^— 7- I — J ^^-^ — j — Vrhms- Three averages of v are shown in Fig. |2j 
Thus, we see that Vrhms is the more representative speed in the ER regime. 
These ideas can be used for a bosonic system as well by replacing t] with ( 
. But C(l) = oo, so Vrhms{boson) = 0, indicative of the liability of bosons to 
condense into a zero momentum state. Although a zero momentum state is 
strictly forbidden for a fermionic system, as t] (1) = In (2), a low momentum 
condensate characterized by the harmonic mean values is possible. The rms 
energy is given by a^E^^^ = y'^^ns + another characteristic value of 

the energy is the velocity averaged momentum a?El =< y"^ >< >= 
ylmsO- + x'^lvlmhs) ■ Obviously, as a ^ 0, E ^ cx). 

Any macroscopic quantity that depends on the momentum distribution 
will be sensitive to the averaging process. This will be even more so in the 
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relativistic regime. In this paper, we look into such effects and investigate 
the possible distribution of the sizes of neutrino structures. In particular, 
we would like to determine whether smaller neutrino structures could have 
formed in the very early Universe. In the next Section, we apply these ideas to 
the gravitating neutrino spheres, and compare the distributions of Keplerian, 
Virial and Jeans and the free streaming scales. 



Neutrino Jeans Mass 



The Jeans mass that is contained within the Jeans radius, Mj = ^^^-^^^ = 
has the momentum dependence jH], 

C (y2 + a;2)7/4' ^ 

7/2 3 

where the constant C = = ^ x 15 x lO^^Mc?), nipi is the Planck 

mass, X = ma = m/T and y = pa = p/T . Any gas cloud of radius greater 
than Rj will contain a mass greater than Mj, in which case gravitation can 
overcome the free streaming motion of the neutrinos to produce a collapse. 
This analysis can be used for any or all the components of the Universe, and 
here we will apply it to the neutrino. These characteristic length and mass 
scales, being dependent on the momentum distribution, are sensitive to the 
averaging process. 

J.R. Bond et al^ calculated Mj by comparing the pressure and the 
density as < Mj >i~ (pf ^)^^^/(p)^5 which was found to peak at Xi = 4.2. 
This expression may be appropriate in the gravitational collapse of a mass of 
a gas against its internal pressure, but the neutrinos are free-streaming since 
decoupling. Thinking that kj is the fundamental quantity that determines 
the Jeans scale, we calculated < Mj >2~ (pf^)^/^/ {kj)^ and found the peak 
to be at Xi = 2.5. We can also find the average of the Jeans mass as < 
Mj >3=< Mj{y) > with peak at 5.0, and represent some characteristic 
values analytically by < Mj >4=< Mj{yrms) > with peak at 4.2, and < 
Mj >5=< Mj{yrhms) > with peak at 1.9. Plot of these values are displayed 
in Fig. [3l The values of x at which these Mj's peak are given in the table 
[TJ The results show quite a wide variations in the Jeans mass ranging across 
a magnitude. Also the time when these structures can form, determined by 
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Table 1: The value Xp at which the different means of neutrino Jeans mass 
peak and also the peak value. As the maximum of the Xp's is 5, the different < 
Mj >'s are caculated at this value and the standard deviation is determined. 



Method 


Xp 


Mj{xp) 


Mj{x = 5) 


Standard Deviation 


< Mj >i 


4.2 


1.75 


1.71 




< Mj >2 


2.1 


0.73 


0.36 




< Mj >3 


5.0 


1.61 


1.61 


0.832 


< Mj >4 


4.2 


2.06 


2.01 




< Mj >5 


1.9 


0.41 


0.21 





X, varies from 1.9 to 5. In essence two groups of neutrino structures could 
appear: the smaller one at x ~ 2 when the neutrino are still quite relativistic 
with Vrhms ~ 0.6, and the ten times larger one at x ~ 5 when the neutrinos 
have become non-relativistic with Vrhms ~ 0.3. 



3 Neutrino structure scales 



To investigate the scales of neutrino structures, let us rewrite the Friedmann 
equation in terms of the Hubble radius = 1/H in the following form, 

K m e \dt) P.0 e e • ^ ^ 

Here the subscript represents the respective values at some reference time 
and pc = is the critical density; p is the total energy density contributed 
by matter and radiation, while 1 — f2 may be considered to be the energy 
density due to curvature, and in particular, /(O = ^Ao^^ + (l~^o)^^ + ^(io^+ 
Q^o + Quresto^/ {^"^ + {v I Xq)) whcrc ^ = x/xq = a/ao = T^/T, ^a^^^^u are 
the respective contributions to the total density at to (in units of critical 
density) by the cosmological constant, pressureless dust, photon, neutrino, 
etc. In the way we have written the neutrino density, we need not make the 
distinction between radiation and matter domination as its density scales 
appropriately in the respective regime 

/ l^uo/e, e^o, (ER) 

The ratio of neutrino and photon density is 
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The temperature ratio T^/T^ is unity at very high temperature and is 
(4/11)^/^ at temperatures below the neutrino decouphng and e"*" — e~ anni- 
hilation, i.e., Tj, < IMeV; the density ratio rises linearly in the NR region. 
The temperature of photon-neutrino equality can be determined by setting 
Eq. ((jl])) equal to unity, whence one finds Xeq = 3.179. 

The most important length scale in cosmological context is the Hubble 
radius; the ratio of the comoving Hubble radii at equality and at arbitrary 
time is given by 



Rn/a J \ Rh J ^ ' 
A length scale relevant to a gravitationally bound gaseous system is the 
radius where the escape speed becomes equal to the average random speed of 
the gas particles viz., v'^ = 2M/R(.sc = ^pRlsc^ ~ smooth 

background of photons will not contribute to the potential well, p = p^. We 
can use the expression for critical density to write the ratio of the comoving 
Hubble and escape radii in a suggestive form that looks similar to Eq. ([5]): 



R. 



where 



fesciO = ^urestoWe + iy/Xoy)/v'. (7) 

We can also define the Keplerian radius Rk by^ = 4^ = op2 with the 
understanding that v stands for the rotational speed; obviously, any particle 
orbiting with Vrot > v^sc will escape. Similarly, the virial radius given by 
2 < v'^ >= M < 1/Rv > leads to = |^ = Another such scale, 

the Jeans radius is given by ^ = = „ o p, . All these length scales are 

some numerical multiples of Resc- 

In any spectral analysis, it is the wavenumber ~ that is fundamen- 
tal, and similarly, in the Jeans analysis, the Jeans wavenumber kj should 
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be averaged over all momenta. We have seen that R^sc is the underlying 
quantity, so let us look into the behaviour of 1/ R^^^. There are a number of 
ways of averaging Eq. (|6]), as discussed in the previous section: 





)> 


^urestO 


Resc / 




e 




)> 


rest 


Resc y 




e 




)> 


^urestO 


Resc / 







2^3/2 '\ , 

from (4) (8a) 



(8b) 

from M (8c) 

where, in the last sub-equation, P stands for pressure and has been 
replaced by 3P/ p. The plot of R^^^c averaged in three ways of Eq. are 
displayed in Fig. [H] for a flat universe, but the general trend for closed and 
open is also the same:i?esc increases to a maximum value, and then decreases 
slowly with the expansion of the universe. But there are drastic differences 
in the values of the last average relative to the first two; e.g., for Q^eq = 

^uea = 0.5 that is shown, we find that the maximum of ( -rk^ ] = 0.757 
occurs at ^ = 0.847, that of (^■^^] = 0.477 occurs at ^ = 0.922, and 



2max 



\ Ir^ ) ~ 1.074 at ^ = 1.632. These different values indicate a spectrum 

V ^ ^0 / 3max 

of possible length scales. We may combine Eqs. |5] and [6] to write 



Rh _ / fesciO /Q^ 
Resc V /^(O 

This relation, displayed in Fig [5l shows that the horizon crossing at 
Rh/ Resc = 1 occurs at different times for different averages. proportional 
to Resci will also show similar behaviour. 

The smallest of all these length scales, -Resc, enters the horizon the earliest, 
while the largest, Rj, enters later. It appears that the time lag between Rj 
entering horizon at the onset of Jeans instabilities, the virialization of the 
cosmic neutrino inside a smaller Ry that entered the horizon earlier, and 
even the Keplerization of the gas in an even smaller length Rk cannot be 
great. 
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We can also calculate the netrino mass within i?esc viz., Me^c = ^puResc 
that simplifies to = ,„ ^^l^^] ,^ , — ttwtj- Here also, as in the previous 
section, the averaging over the momentum of M^sc can be done in many 
difi^erent ways. Below, we give three which we think are useful: 



Mesc\ 0.5 / {y/xofe 



'Mesc\ 0.5 / {y/x^f \ I {y/xQ)( 



1/2 



(10) 



(11) 



M,,,\ _ 0.5 /e + {y/xofy / {y/xo)e V^' ^^2) 



\Rho/s yAh^ZroX {y/xo? I \(e + (2//xo)2)3/^ 

A plot of {M(.sc/ Rhq) is shown in Fig. H] for three different averages. 
Although the nature of the curves are same, the peaks are seen to occur at 
different positions. 



4 Virial Equation and Moments 

Moments of the virialized quantities have played very important role in re- 
vealing the features of large astrophysical objects. Some of these methods 
are used here to describe the large scale neutrino structures that should be 
the precursors of the galactic and the supercluster scales. The formation of 
large-scale structures proceeds with the gravitational collapse of clouds of 
matter composed of baryons, cold dark matter and hot dark matter whose 
major constituent is the neutrino. Much work has been done to investigate 
the evolution of size and eccentricity of spherical and spheroidal / ellipsoidal 
clusters [Hl[9l[l0], and and also the relavant masses. In the majority of these 
works, Virial theorems / equations and their moments of various orders have 
been used. In this section, various moments of the Jeans mass are calculated 
using using the ideas of Virial moments. The results are compared with the 
values presented in the previous section. 

In a system of N particles, gravitational forces tend to pull the system 
together and the stellar velocities tend to make it fly apart. It is possible 
to relate kinetic and potential energy of a system through the change of its 
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moment of inertia. In a steady-state system, these tendencies are balanced, 
which is expressed quantitatively through the Virial Theorem. A system 
that is not in balance will tend to move towards its virialized state. The 
Scalar Virial Theorem tells us that the average kinetic and potential energy 
must be in balance. The tensor Virial Theorem tells us that the kinetic 
and potential energy must be in balance in each separate direction. The 
scalar Virial Theorem is useful for estimating global average properties, such 
as total mass, escape velocity and relaxation time, while the tensor Virial 
Theorem is useful for relating shapes of systems to their kinematics, e.g. the 
flatness of elliptical galaxies to their rotational speed. 

The Virial Equations of the various orders are, in fact, no more than 
the moments of the relevant hydro dynamical equations. The scalar Virial 
Equation for a system is given by 

1 d^I 

~-^I<^U (13) 

where the moment of inertia about the origin I = ^ f pr'^dV, the kinetic 
energy is = Jypv'^dV and the potential energy is U = ^ Sv i 
where pcis the core density. Thus the Virial Equation of the 1st order is 
given by 

V V V 

The Virial Equation of the 2nd order is given by just multiplying the 
integrand by r before integrating: 



1 [ , I (p 

V V V 

Similarly, the higher order equation may be written. For the steady state, 
2K + U = 0. This gives the virial radius of the spherical cluster: 

Writing p and v in terms of x = ma = m/T and y = pa = p/T, as 
discussed in Introduction section, we have 



9 



Rv = Constant. — ^ ( — ^ (17) 

and hence the expectation value of this virial radius is given by 

oo 

/y^ / X \ ^ 

dy =77 ( — ) (18) 

Its evolution has been shown in [71 

The Virial Mass of the cluster contained within this radius may be cal- 
culated as < My >= constant. < Ry >^< p > . Its variation with x is as 
expected ([9]), but peaks at x = 3.73 in contrast to the position of peaks of 
Jeans masses calculated in section [2l where the peaks have occurred at x = 
1.9, 2.1, 4.2 and 5.0. [HIT] 



4.1 Virial moments of Jeans mass 

As stated earlier, in the virial method, we take the moments of the equation 
of motion. These equations obviously involve the moments of the distribution 
of density, pressure, velocity, gravitational potential, etc. Here we are taking 
the spatial moments of various orders of Jeans mass given by [2j The 1st 
moment of Jeans mass is given by 



oo 

f 

(MjRv) = constant. / dy- —{MjRi 





oo 



fi 4 

yy X 
dy 777 (19) 
{ey + 1) {y^ + x^f^ 

The 2nd and 3rd moments may be written in similar fashion. The vari- 
ations of the 1st, 2nd and 3rd moments with x have been shown in [SlThe 
1st and the 2nd are found to peak and at x = 12.2 and 27.02 respectively, 
but the 3rd moment appears to plateau off from x ~ 100 to large values of 
X. Jeans masses may be calculated by dividing these moments by < Ry >, 
< Ry > and < Ry > respectively. The variation of these masses with x 
has been shown in M (For comparison < Mj > and < My > have also been 
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Table 2: Temperatures at which the peaks of the Jeans masses occur for 
different neutrino masses. 



Mj from 


X for Mj peaks 


Temp. (K) corresponding to x for 


m,y= 1 eV 


m^= 0.2 eV 


m^= 0.01 eV 




1.9 


6100 


1220 


61 


\p) 

(Kjf 


2.1 


5519 


1104 


55 


< Ry >^< P > 


3.73 


3107 


621 


31 


Mj{yrrns) 


4.2 


2760 


552 


28 


{Mj{y)) 


5 


2318 


464 


23 


1st moment of Mj 


6.48 


1789 


358 


18 


2nd moment of Mj 


8.01 


1447 


289 


14 


3rd moment of Mj 


9.58 


1210 


242 


12 



plotted). They have peaked at x = 6.48, 8.01 and 9.58. Following table gives 
the corresponding temperatures at which the peaks of the Jeans masses occur 
for a few typical neutrino masses: 

From the above analysis, it is seen that the large-scale structures of neutri- 
nos of different mass and random velocity distribution can form at different 
neutrino temperatures, corresponding to different time. The earliest peak 
that occurs at x = 1.9 corresponds to the time when the neutrino temper- 
ature was 1220 K for the 0.2 eV neutrino. Similarly, the latest peak occurs 
at X = 9.58 corresponding to a temperature of 242 K. In between these two 
values, it is found that the Jeans mass peaks at a number of different x. Thus 
it can be interpreted to mean that a distribution of neutrino structures of 
different masses and of different ages should be in existence. Typical masses 
of these structures range from 6 x 10^^ to 4.5 x IO'^^Mq. 

5 Free Streaming 

It is well known that free streaming of collisionless particles wipe out any 
structures that form. In this context, there are other two scales that are 
relevant. The first is the particle horizon Rp occurring at the distance to 
which light would travel since the big bang. The comoving particle horizon 
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can be written as 

Rp di' 



, (20) 

No massive particle can cover a distance larger than this. The other is the 
freestreaming length Rp^ the distance a particle travels since the big bang; 
it is given by an expression similar to Eq. (1201) . but including the particle 
speed V inside the integral: 



Jy/x 



-z I ^ rest eq ^ , ^ '7 eg ^'■u rest eq 



Any scale that is smaller than Rp will be wiped out by free streaming. 
The integral Eq. (1211) in our case, for different values of fi^eg, can be written 
as Jacobian elliptic functions |I6] 



Rf 



= CF{(j)/k) (22) 

flat 



Considering the flat case where = 0.5, the solution is = ^(y/^^i) — 

76? ' -, I 3C(3) y T^Y^ 

sin^ = and = ^^^-^^ ■ A plot of (/?i./ei?/.. J//at is shown 

J^iq+y^+y y/'-^iq+y^+^eq 

in Fig. - [10] for three different mean values of y. This free streaming length 
is seen to saturate to the values of 4.193, 3.899 and 2.555 as x — )■ 00 , when 
evaluated with yrms, Vmean and yrhms respectively. 

For non-fiat universe, Rp/^Rn^^ is given by Eqs. f l2T]) and f l22l) with = 

M^Jl + {y/x,,)\ tan2(0/2) = 2P = 1 + 

(2//xeg)2 (i^) v^l + {y/xeqf. Plots of Rp/^Rh.^ for fi^eg = 0.45, are 
found to saturate to 2.883, 2.113 and 1.690, for the three means of y; and 
in the closed universe of Vt^eq = 0.55, the free streaming length saturates 
to 2.683, 1.908 and 1.495 for the respective means of y. The nature of the 
curves are similar to that of fiat ones. 



6 Discussion and conclusion 



12 



From above analysis, it is seen that the large-scale structures of neutrinos of 
different mass and random velocity distribution can form at different neutrino 
temperatures, corresponding to different time. The earliest peak that occurs 
at X = 1.9 corresponds to the time when the neutrino temperature was 1220 
K for the 0.2 eV neutrino. Similarly, the latest peak occurs at x = 9.58 
corresponding to a temperature of 242 K. In between these two values, it 
is found that the Jeans mass peaks at a number of different x and hence 
it can be interpreted to mean that a distribution of neutrino structures of 
different masses and of different ages should be in existence. Typical masses 
of these structures range from 6 x 10^^ to 4.5 x 10^*^ Mq. These are in the 
order of very large super-clusters which formed first in the hot dark matter 
scenario. But we should note that as the Universe cools down, and x becomes 
large, the scale of these structures will come down because the leading term 
in the Jeans mass, Eq. ([2]) goes as x~^^'^. Thus for 1 eV neutrino at the 
present temperature of 1.96 K (~ 10~^e\^), the Jeans mass comes down to 
the galactic scale. Accordingly, the large superclusters that formed initially 
would have fragmented to the galactic size. 

Thus the massive neutrino can play a very important role in the forma- 
tion of structures. These particles decouple from the rest of matter at a very 
early time when the temperature was 1 MeV. Since then they have been 
cooling down independently without interacting with other matter except 
gravitation. So the evolution of this component can be considered to occur 
separately. We have done the Jeans analysis taking into account only neu- 
trino. The results show the possibility of different sizes of neutrino structures 
forming at different time. In particular, two distinct sizes make their appear- 
ance: the smaller ones which could form at relatively higher temperature of 
T ~ m/2 when the neutrinos were still relativistic, and the ten time bigger 
ones that form when T ~ m/5. As these neutrino structures have a mass 
of order of m^^/m^, these are indeed very large, in the scale of large super- 
clusters. The Jeans mass at lower temperature is of the order of galaxy. In 
this work we have considered the neutrino as a totally independent com- 
ponent of the Universe that interacts with nothing except gravitation with 
each other. In reality, these neutrino interact gravitationally with the rest of 
the matter. So in any volume, there is more matter than just neutrino that 
is trying to bring about a gravitational collapse. Hence the actual scale at 
which the neutrino could begin to collapse should be smaller than the value 
we have derived. 
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Figure 1: Momentum distribution of the free cosmic neutrino as a function 
of y = pa = p/Ty. The distribution is maximum at Umax 3.131, and other 
characteristic values are Urms =< Z/^ >^/^= 3.597 and Urhms =< 1/?/^ >~^/2= 
1.613. Such momentum distribution give rise to a spectrum of neutrino 
structure scales. 
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2 4 6 8 10 

Figure 2: The average speed of the massive neutrino as a function of x. The 
sohd curve is Vmean-, the dashed curve is Vrms and the dotted curve is Vrhms- 
In the extreme relativistic regime as x — > 0, it is seen that Vmean Vrhms] 
towards the non- relativistic end it is closer to tVm.s- 
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<Mj> 



2 . 5 




Figure 3: The distribution of Jeans Mass. The dashed curve represents the 
value calculated by Bond, et al, i.e., < Mj >i The dotted curve is < Mj >2, 
the thick solid curve is < Mj >3, moderately thick curve is < Mj >4 and the 
normal curve is < Mj >5. Two distinct groups of Jeans mass are seen, one 
that peaks at x ~ 4 to a value of ~ 4, and another that peaks at x ~ 2 to 
little less than 1. The first group appears to be related to the rms speed and 
the second to the rhms speed. As pointed out in the text, at the relativistic 
regime it is the rhms speed that is more representative. So the smaller second 
group of structures can start to form with smaller mass at earlier times. 
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Figure 4: The distribution of Mass contained within escape radius calculated 
by different averaging process given in equations [101 [ID and [121 The dashed 
curve represents the value calculated by the 1st equation, the dotted is that 
of 2nd equation and normal curve is from the third one. 




Figure 5: The ratio of Rh/Rssc for the three averages. The third average 
enters the horizon the latest. 
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Figure 6: The ratio of Resc/^Rn for the three averages. 




Figure 7: Evolution of Virial Radius < Ry > of a cluster. 
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Figure 8: Virialized moments of the first three orders of Jeans masses. The 
first moment peaks at x ~ / T = 6.48. The second moment peaks at x = 
8.01 and the 3rd one at 9. 58. The sohd, dashed and dotted hnes are for 1st, 
2nd and 3rd order moments respectively. They are appropriately scaled to 
make of equal heights for the comparision purpose. 
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Figure 9: The average cluster masses (Jeans and Virial) against x. The 
normal curve is for virial mass, thick continuous curve for Jeans mass, the 
normal dashed curve for the Jeans mass calculated from the 1st moment of 
the mass and normal dotted and thick dotted are those from 2nd and 3rd 
moments respectively. 
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Figure 10: The ratio of Rf/S,Rho for the three averages for the flat case. 
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